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Abstract 

We investigate the i/e — > oscillation in the framework of three generations when neutrinos pass 
through the earth. The oscillation probability is represented by the form, P(i^e i^n) = A cos 5 + 
B sin 5 + C in arbitrary matter profile by using the leptonic CP phase S. We compare our approximate 
formula in the previous paper with the formula which includes second order terms of a = Am^i/Am^i 
and Si3 = sin 013. Non-perturbative effects of a and sia can be taken into account in our formula and 
the precision of the formula is rather improved around the MSW resonance region. Furthermore, we 
compare the earth matter effect of A and B with that of C studied by other authors. We show that 
the magnitude of A and B can reach a few ten % of C around the main three peaks of C in the region 
E > 1 GeV by numerical calculation. We give the qualitative understanding of this result by using our 
approximate formula. The mantle-core effect, which is diff'erent from the usual MSW effect, appears 
not only in C but also in A and B, although the effect is weakened. 

1 Introduction 

The first evidence of neutrino oscillation have been discovered in the atmospheric neutrino experiments 
and the mass squared difference jAml^l and the 2-3 mixing angle 623 [J have been measured. Also the 
deficit of solar neutrino strongly suggests the neutrino oscillation with the Large Mixing Angle (LMA) 
solution for Arn^i and 9i2 [2]' This has been confirmed by the KamLAND experiment by using the 
artificial neutrino beam emitted from several reactors On the other hand, only the upper bound 
sin^ 2^13 < 0.1 is obtained for the 1-3 mixing angle @|. Thus, the values of the mass differences and the 
mixing angles are gradually clarified. Our aim in the future is to determine the unknown parameters like 
the sign of Am|i, 813 and the leptonic CP phase 6. 

The simple analytic formula for estimating the matter effects is useful in order to study these pa- 
rameters because neutrinos pass through the earth in most experiments and receive the matter potential 
represented by a = V2Gi?iVe, where Ne is the electron number density and Gp is the Fermi constant. In 
the case of short baseline length, we can approximate the density as constant because the variation of TVg 
is small. However, the longer the baseline is, the larger the matter effect is. In previous papers, several 
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approximate formulas have been proposed in order to include the effect of varying density. Classified by 
the neutrino energy E, there are following approximate formulas: low energy formulas by the expansion 
in the small parameter 2Ea/\I^m\i\ <C 1 or S13 = sin^ia <C 1 high energy formulas by the expansion 
in /^m\^/2Ea ^ 1 or a = /^m\^//S.m\^ <C 1 6 , and the formulas by the expansion in 2E5a/ /^m\i ^ 1 
[7], where 5a is the deviation from the average matter potential. 

On the other hand, there is the method to approximate the earth matter density as three constant 
layers in the case of mantle-core-mantle 0. It was discussed in refs. ^ how the probability is enhanced 
when neutrinos pass through periodically varying density. Then, it was pointed out in ref. |10| that the 
mantle-core effect, which is different from the usual Mikheyev-Smirnov-Wolfenstein (MSW) effect 
appears in the oscillation probability. More detailed analysis has been in refs. |12l This effect is 
interesting because the large enhancement of the probability can occur even if both the effective mixing 
angles in the mantle and the core are small. In recent papers the possibility for measuring the ^13 
in atmospheric neutrino experiments has been discussed by using this mantle-core effect. They concluded 
that the value of ^13 can be measured in some cases. 

In our previous papers, we have shown that the oscillation probability for v^. v^^ transition is repre- 
sented by the following form, 

P{ve Vfj^) = A cos (5 + B sin (5 + C (1) 

in constant matter jj^ and also in arbitrary matter By using this general feature for the CP depen- 
dence, the method for solving the parameter ambiguity problem pointed out in refs. |17| is discussed in 
ref. JS]. Each coefficients has an order A = 0(si3a), B = 0(si3a) and C = 0(5^3) + 0{a^) on the two 
small parameters a = Am|]^/Am§]^ ~ 0.04 and S13 ~ sin 6*13 < 0.2. In the case of a < S13, the ratio of 
A, B to C are given by A/C — O^a/sis) and B/C — O^a/sis). So, it is expected that the CP violating 
effect due to A and B becomes large and can reach a few ten % of C even for the case that neutrinos 
pass through the earth core. However, the effect due to the CP phase has not been taken into account in 
previous works. 

In this paper, as the preparation of studying earth matter effect, we review our approximate formula 
introduced in ref. ^Hl 

A^ 2c23S23Re[Sf:,S'^,], (2) 
B c ~2c23S23MS%S!^,], (3) 
\Slfcl, + \S!^fsl„ (4) 

where and S!^^ are the amplitudes calculated from two-generation Hamiltonians H( and Hh- He is 
represented by Am^i and 612 and Hh is represented by Am^^ and 613. We show that our formula includes 
the non-perturbative effect of a and S13 and the precision is rather improved around the MSW resonance 
regions compared to the well known simple formula |20[ 12 1| . which includes up to second order terms of 
a and S13. Furthermore, we compare the earth matter effect of A and B with that of C by using the 
Preliminary Reference Earth Model (PREM) |22| in the case of two reference baseline length. We show 
that the magnitude of A and B can reach a few ten % of C around the main three peaks of C in the region 
E > I GeV by numerical calculation. This means that the above perturbative estimation is valid even in 
the case of including non-perturbative effect. We give the qualitative understanding of this result by using 
our approximate formula. The mantle-core effect, which is different from the usual MSW effect, appears 
not only in C but also in A and B, although the effect is weakened. 
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2 General Formulation for Neutrino Oscillation Probabilities 



In this section, we review the exact formulation of neutrino oscihation in arbitrary matter profile based 
on ref. At first, let us parametrize the Maki-Nakagawa-Sakata (MNS) matrix U 23 , which connects 
the flavor eigenstate Va with the mass eigenstate Vi, by the standard parametrization |24j 

U= 023rsOi3TlOi2, (5) 

where Ts = diag(1, 1, e'*) and Oij is the rotation matrix between i and j generation. As the matter 
potential only appears in the (ee) component of the Hamiltonian, O23 and Tg can be factored out. So, we 
can rewrite the Hamiltonian in matter as iJ = O23T s H' sY ■ H' is the reduced Hamiltonian defined 
on the basis v' = (023ri5)^i^. H' is a real symmetric matrix and the concrete expression is given by 

H' = Oi3Oi2diag(0,A2i,A3i)(Oi3Oi2)^ + diag(a(i),0,0), (6) 

where Ay = /\ra1^/2E. The number of parameters in the Hamiltonian H' is fewer than that in the 
original Hamiltonian H by two. This is useful to calculate the oscillation probability simply. If we define 
the amplitude S'^^ of v'^ v'^ transition as {a(3) component of time ordered product 

(7) 

the oscillation probability is given by 

A cos J + S sin (5 + C, (8) 

2c23S23MS'leS're\. (9) 
-2c23S23lHS'leS'rel (10) 
l^;eP4 + l^;el'4 (11) 

as in 16' . /^From the eqns. l(^- Hll|) . one can see that the probability for v^, transition is represented 

by two components of the reduced amplitude, S'^^ and S'^^. Namely, the matter effect for the oscillation 
probability is only contained in the two components. 



S' ^ Texp 



H'{t)dt 



A = 
B ^ 
C = 



3 Non-Perturbative Effect in Our Approximate Formula 

In this section, we numerically calculate the amplitudes S'^^ and S'^^ introduced in the previous section 
by using the PREM. Then, it is explained how we obtain the hint for the basic concept on deriving our 
approximate formula. As an example, the approximate formula in constant matter is derived explicitly and 
is compared with the formula in refs [211 1^ , which includes up to second order of the small parameters 
a = A21/A31 and S13. As a result, it is shown that our approximate formula includes the non-perturbative 
effect which becomes important around the MSW resonance. 

3.1 Behavior of Reduced Amplitudes in Earth Matter 

Let us calculate the amplitudes S'^^ and S'^^ for the case that neutrinos pass through the earth. We 
use the PREM as the earth density model and we choose two reference baselines, 6000 km and 12000 
km. Fig. 1 shows how the matter density changes along the path of neutrinos. In fig. 2, we plot the 
values of the amplitudes S'^^ and S'^^ corresponding to the neutrino energy 0.03-20 GeV. Here, we use the 
parameters Am^^ = 7 x 10^''^eV^ and sin^ 2^12 — 0.8 as indicated from the solar neutrino experiments 
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and the KaniLAND experiment, Arn^i = 2 x 10 ^eV^ from the atmospheric neutrino experiments and the 
K2K experiment, and sin^ 20i3 = 0.1 within the upper hmit of the CHOOZ experiment. 
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Fig. 1. Matter density in the PREM with basehne length 6000 km and 12000 km from left to right. 
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Fig. 2. Energy dependence of S'^,, and S'^^ with baseline length 6000 km and 12000 km by using the 
PREM. The solid and dashed lines represent S'^^, and S'^^ respectively. 

It is found from fig. 2 that S'^^ and S'^^ become large in low energy and high energy, respectively, for 
both baselines and the regions, where S'^^ and S'^^ dominantly contribute, are separated to each other. 
In other words, the MSW effect related to the 1-2 mixing and 1-3 mixing angles are mainly included in 
S'^g and S'^g, respectively. We have derived the approximate formula for arbitrary matter profile by using 
this feature in ref. 19 . Concretely, and S'^^ are calculated from two kinds of different Hamiltonians, 
which are given by the 1-2 and 1-3 subsystems, respectively. 



3.2 Procedure of Deriving Approximate Formula 

The idea introduced in the previous subsection is actually realized as follows. We use the two small 
parameters a = A21/A31 and S13. Then, our approximate formula is calculated by the following three 
steps. 
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1. We define two Hamiltonians in the 1-2 and 1-3 subsystems taking the limit of S13 — > and a ^ in 
® as 



A2isl2+a{t) A21C12S12 

A21C12S12 A2lC?2 

A31 



Hh = 



/ A3is?3 + a(i) A31C13S13 



y A31C13S13 AsiC^g 



(12) 



2. We calculate two amplitudes and S'' from the Hamiltonians iJ^ and by the equations 



S-^ = Texp 



-i / Hi{t)dt 
Jq 



, 5'''= Texp 



Jo 



(13) 



3. We replace the amplitudes in (O-dJ as S'^^ Sf^^ and S'^^ S'^g. 

3.3 Approximate Formula in Constant Matter 

Next, let us review the approximate formula in constant matter based on ref. |19j . According to the 
procedure in the previous subsection, we substitute the Hamiltonian Hg in constant matter given by p2ll 
into (|13() and we obtain S^^ as 



S'^g = [exp{-iHiL)] = -J sin 26'f sin exp 



. A21 -f- a 



(14) 



where cjje = A(L/2 and the subscript £ represents the quantities calculated from Hi. The concrete expres- 
sions for the mass squared difference and the 1-2 mixing angle in matter are given by 



Af sin 2^12 
A21 sm26( 



cos 26*12 - -r — I + sin^ 26*12. 
A21 , 



(15) 



These are well known expressions in the framework of two generations. The contribution of the low energy 
MSW effect, which is dominant around the energy region determined by a ^ A21 cos 26*12, is included 
in mainly 5^^. The phase factor in H14() does not contribute when we calculate the probability in two 
generations. However, this gives important contribution on the calculation of the terms dependent on the 
CP phase in three generations. 

Similarly, we obtain S!^^ by substituting Hh in constant matter given by H12|) into (|13|l as 



S'^g = [exp(-ii//iL)]^g = -i ain20h sin (l)hexp ^-i— ^y!^L^ , 



(16) 



where (j>h = A/iL/2 and the subscript h represents the quantities calculated from Hh. The concrete 
expressions are given by 



Ah sin 26* 



13 



A31 sin 26*^ 



cos 26*13 - ) +sin^ 26*13. 
A31, 



(17) 
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One can see that these expressions correspond to those obtained by the replacement A21 A31 and 
6*12 Oi3 in (|14|l and (|15|l . The contribution of high energy MSW effect, which is dominant around the 
energy region determined by a '--^ A31 cos 2^13, is included in mainly S!^^. We can calculate A, B and C in 
constant matter as 

P{ve-^ Vi,) = Acoad + Ba\n5 + C, (18) 

A32i 

A ~ sin 20f sin 2023 sin sin sin cos — ^ — i (19) 

A32L 

B ~ sin 26'f sin 2^23 sin 20/1 sin 0£ sin </);i sin — ^ — , (20) 

C ~ C23 sin^ 20f sin^ 0£ + S23 sin^ 20/1 sin^ (21) 
from these expressions. These approximate formulas are similar to the following well known formulas 



A21A31 aL (a-A^AL A31L , ^ 

A ~ — — -Ci3sm26'i2sm26'23sm26'i3sm — sm cos — — , (22) 

a(a — A31 ) 2 2 2 



A21A31 . aL . (a-A3i)L . A31L 

B ~ — — -ci3 sm 26112 sm26'23sm26'i3sm — sm ^ sm— — , (23) 

a[a — A31) 2 2 2 

C 4^4sin2 20i2sin2^ + ^^^4sin^20i3sin2^^^^. (24) 

2 (a — A31 )^ 2 



These formulas are often used in order to analyze the property of neutrino oscillation because they have 
very simple form and approximate the exact values with a good precision. 

In the following, we compare the probability calculated from our approximate formula (|19|I - H21() with 
that from the formula H22() - (|24|l in the case of constant matter. We calculate P(i^e — > ^^/i) for two kinds 
of baselines, 3000 km and 6000 km. We use the parameters sin^ 26*23 = 1 and 5 — 0° in addition to those 
introduced in fig. 2. Furthermore, p = 4.7g/cm'^ and Vg = 0.494 are used as the matter density and the 
electron fraction. The result is given in fig. 3. 




Fig. 3. Comparison between the probabilities P{ve — >■ i^ii) calculated from our approximate formula 
(Cnil-lEU and the formula |221)-(I2i|l- Two baseline length are chosen as 3000 km and 6000 km from left to 
right. The solid, dashed and dotted lines are the exact, our approximate formula and that from H22|) - (|24|l . 



6 



It is found that our approximate formula has good coincidence to the exact one even around the MSW 
resonance, compared with that from the formula H22|l - H24|) . We consider the reason for the difference in 
the next section. 



3.4 Comparison of Approximate Formulas 

Let us explain the order counting of a and S13 in our approximate formula. In the limit a — > 0, we 
obtain S'^^ = 0. Therefore, we can write down the order of S'^^ as 



S" — 



0{a) + 0{a^) + Oia^) + ■ ■ ■ 
+ Oisua) + 0{si3a'^) + 0{si3a^) + ■ ■ ■ 
+ 0(sl,a) + Oisl^a') + Oisl^a') + ■ ■ ■ 

+ • (25) 

Note that a is included in all terms. Here, if we take the S13 — > 0, only the first line is remaining. In these 
terms, all orders of a are included and the first line considered to have a larger contribution compared to 
the following lines, because of the increasing exponent of S13. This is confirmed by the comparison with 
the exact formula in fig. 3. In the same way, we obtain = in the limit S13 —^ 0. Therefore, we can 
write down the order of S' as 



0{s,3) + 0{sl,) + 0{sl,) + -- 
0(si3a) + 0(s?3a) + 0(s?3a) 



+ Oisiaa^) + 0{sl^a^) + Oisl-^a^) + ■ ■ ■ 

+ • (26) 

Here, if we take the limit a — > 0, only the first line is remaining. In these terms, all orders of S13 are 
included and the first line is considered to have a larger contribution compared to the following lines, 
because of the increasing exponent of a. 

Our method includes both, the terms of higher order of a in (|25|l and also those of S13 in H26|l . So, 
this new approach is not a systematic expansion. However, our method is not in contradiction to the well 
known formula l|22|l - l|24|) . which takes only the first order term of a and S13 in H25|l - 126|) regarding them as 
small parameters. In addition, higher order terms of the perturbative expansion, which are not included 
in the formula (|22() - (|24|l . are now also included in our formula. 

In the following, let us investigate the difference between these two methods more concretely. As seen 
in fig. 2, the contribution of S'^^ is dominant in the energy region E > 1 GeV. So, we can roughly consider 
as 



P{i^e ^ i^fj.) — C ~ S23 sin^ 26h sin^ 



(27) 



Note that we have the relation (|17|l between the mass squared differences and the mixing angles in vacuum 
and in matter. Expanding the right hand side of (|17f) on the mixing angle in vacuum, we obtain 



Ah_ _ sin 26*13 
A31 sin 29h 



a 



2aA 



31 



(A31 - a) 



2 ''13 



2(A3i - a) 



4 ■'13 



(28) 



The condition for convergence is given by 



4aA3isf3 
(A31 - a)2 



< 1. 



(29) 
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This condition is not satisfied around the MSW resonance region, where A31 ~ a. Namely, the pertur- 
bative expansion becomes inconvergent. However, substituting the above expression H28|l into H27|l of the 
oscillation probability and taking only the first term, we obtain 



'29 



23" 



(A, 



13 . 2 

sm ■ 



^31 



31 



(30) 



It gives the finite value in the limit A31 — > a as 

P(zye^l.^)~4c?3(si3A3lL)^ (31) 

This is due to the product of the infinity of effective mixing angle and the zero of the effective mass squared 
difference in the probability (|30|l . If we take the limit A31 — s- a directly in the non-perturbative expression 
(|27|l . we obtain 

P{Ve V^) = s23C?3 sin2(si3A3iL). (32) 

We find that the difference between the perturbative formula (|31|l and our formula l|32|) is the sin factor. 
In the case of the short baseline length L, the perturbation gives a good approximation, but the longer 
the baseline L is, the worse the perturbation becomes, as shown in fig. 3, although the probability has a 
finite value. Concretely, if the condition 

L < (33) 
S13A31 

is satisfied, the perturbation gives a good approximation. Around the MSW resonance region, the pertur- 
bation breaks down because the coefficients of the higher order terms a or S13 become large. Therefore, 
it is needed to involve the higher order terms of a and S13, in order to make a good approximate formula. 
Our method partially realizes this request. 

At the end of this section, let us give a brief comment. In ref. |21)[I2H . the formula calculated by single 
expansion on a is also given and this includes all order terms of S13. So, this approximate formula gives 
a good approximation in the high energy MSW resonance region compared with the formula 122(1 - (|24|l . 
while the difference between the single expansion formula and numerical calculation becomes large in the 
low energy region as commented also in ref. |21| . 

4 Earth Matter Effect for A, B and C 

In this section, we perform numerical calculations of A, B and C by using the PREM. We give a 
qualitative understanding of the behavior of the coefficients A, B and C by matter effects of the mantle 
and the core. 

4.1 Numerical Calculation of A, B and C 

In the previous section, the order of the reduced amplitudes are estimated as S'^^ = 0{a) and S'^^ — 
0(513) in the case that we takes only the first order term of a and S13 in H25|l - H2()l) . The order of coefficients 
are also obtained as A = ©(sisa), B — 0(si3Q!) and C = 0{s\^) + 0[a'^) by substituting S'^^ = 0{a) and 
S'^g = 0{si_z) into ((TT|l - lfT^ . In the case of a < S13, the magnitude of ratios is give by A/C = 0(q;/si3) 
and B/C = 0(a/si3). Therefore, it is expected from the perturbative point of view that the CP violating 
effect due to A and B becomes large and can reach a few ten % of C. However, because non-perturbative 



8 



effect becomes important in MSW region as shown in fig. 3, the CP violating effect should be investigated 
more carefully. 

At first, let us numerically calculate how the coefficients A, B and C are enhanced by the earth 
matter effect, in the case of the baseline length L = 6000km and 12000km for sin^ 26'i3 = 0.10, and 0.04, 
respectively. These values of sin^ 2^13 correspond to the values within the upper bound of the CHOOZ 
experiments. The PREM is used as the earth matter density and the same mass squared differences and 
the mixing angles given in sec. 3 are also used. 
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Fig. 4. Energy dependence of A, B and C by numerical calculation. We use the PREM as earth matter 
density model with two baseline length 6000 km and 12000 km, and we choose the 1-3 mixing angle 
sin^ 26*13 = 0.10 and 0.04 as representative values. 

In the case of L = 6000 km, the behavior can be understood by using the formulation H19|l - H21|) in 
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constant matter. The value of C becomes large around E = 5 GeV, which comes from the enhancement 
of the effective mixing angle sin Oh for a ~ Ah, and then oscillates depending on the factor sin (j>h. The 
values of A and B become small compared with that of C in high energy region, as the suppression factor 
S'^g (X 1/E. See details in refs. |15[ I2U[I^ for example of the constant matter density. 

In the case of L = 12000 km, three main peaks appear in C. On the other hand, the pattern of the 
enhancement for A and B seems to become more complicated than that of L = 6000 km. In the next 
subsection, we give a qualitative understanding of the above results by using our approximate formula for 
A, B and C. 

We also represent the values of A, B and C around the energy of the three peaks of C in Table 1. These 
values are computed by the numerical calculation using the PREM. 



L (km) 


sin^ 2^13 


E (GcV) 


A i\A/C\) 


B {\B/C\) 


C 


peak type 


6000 


0.10 


4.9 


-0.025 


(8.3%) 


-0.003 (1.0%) 


0.301 


MSW(mantle) 


6000 


0.04 


4.8 


-0.017 


(12.0%) 


-0.001 (0.7%) 


0.142 


MSW(mantle) 


12000 


0.10 


2.1 


-0.061 


(17.9%) 


0.066 (19.4%) 


0.340 


MSW(core) 


12000 


0.10 


2.8 


0.017 


(6.3%) 


-0.028 (10.4%) 


0.269 


mantle-core 


12000 


0.10 


5.2 


-0.026 


(7.4%) 


0.013 (3.7%) 


0.352 


MSW(mantle) 


12000 


0.04 


2.0 


-0.068 


(31.6%) 


0.038 (17.7%) 


0.215 


MSW(core) 


12000 


0.04 


3.0 


-0.030 


(6.0%) 


-0.037 (7.4%) 


0.500 


mantle-core 


12000 


0.04 


5.4 


-0.014 


(8.9%) 


0.015 (9.6%) 


0.157 


MSW(mantle) 



Table 1. Resonance values of A, B and C calculated numerically by using the PREM with baseline 
length 6000 km and 12000 km, sin^ 29i3 = 0.10 and 0.04. 

Table 1 shows that the coefficients A and B can be rather large at the three main peaks of C. The 
absolute values of A and B become about 0.06 at the peak of the core. Furthermore, the ratios |^/C| and 
\B/C\ also become a few ten % even for the case of including non-perturbative effect. These energy regions 
E — 2 6 GeV are explored by the atmospheric neutrino and the long baseline experiments. In actual 
experiments averaging of various parameters are necessary for example, energy, zenith-angle distribution, 
the sum of particle and anti-particle, and so on. Therefore, the CP phase effect may be weakened to some 
extent, but we consider that the CP phase effect should be estimated precisely in order to determine the 
value of 013 in future experiments. 

4.2 Approximate Formula in Matter with Three Layers 

In order to give a qualitative understanding of the results obtained in the previous subsection, let us 
approximate the earth matter density with baseline L = 12000 km as three constant layers such that the 
first and the third layers have the same density and length. 

At first, we calculate the amplitude from the low energy Hamiltonian He- We use the superscript 
m and c for representing the amplitude in the first and third layer (mantle), and the second layer (core). 
Taking the limit S13 0, the amplitudes 5*™, S'^ and so on vanish. Only four terms contribute to the 
amplitude in three layers Sf^^ as 

n£ c»m nc nm , nm nc c»m j c»7n nc nrn , nm qc om ('^A\ 

fie ^ fie'^ee'^ee ' fifi'^ fie^ ee ^ fie'^ efi^ fie fifi^ fifi^ fie' V*^ / 

Substituting ((TUl and 

= [exp(-ifff i)]^^ = (cos 4>T + i cos 20^ sin C) exp f J , (35) 
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SZ = [eM-^HrL)l, = (cos C - I cos 20™ sin C) exp ( i ^ L"' ) (36) 



into we obtain 

^^e = --XP (-^ ^"^ + '°7"'"'"^'^0 Fic^T,'^r,07,0d, (37) 

where the function F is defined by 

^(0™, 0"\ 0") = sin 20™ cos 0^= sin 26*'" + cos^ 0" sin c?!)'^ sin 261'= + sin^ 0" sin cf)" sin(26l= - 461'"). (38) 

We can easily extract the physical meaning from this expression, although this function F becomes the 
same one as given in refs. ^| after a short calculation. We describe the meaning of each term later. 

Next, we calculate the amplitude Slf,, taking the limit a 0. In this limit, S'™ , SJ^^ and so on vanish, 
so the amplitude S!^^ in three layer is calculated as 

nh nm oc om i om oc nm , nm nc nm i nm nc nm f '^Q\ 

^re '-'re'-'ee'-'ee ' '-'rr^Te^ee ^Te^er^re ^ ^tt'^^tt^tb' \'^^ ) 

Substituting (|16|l and 

SZ = [eM-iHJ^L)le = (cos C + i cos 2^ sin C) exp -z^i^L™ , (40) 



5™ = [eM-^HJ^L"^)]rr = (cos C - * COS 20™ sin C) exp [~t ^ L^j (41) 
into we obtain 

^:'e - --XP F(c,<^^;^r,^D, (42) 



which corresponds to equation (|37|l by replacing the subscript and superscript as {£) —f (h). 

Substituting H37|l and (|42|l into (|5))- (|ll|l . the coefficients A, B and C in three layers are given by 

A ^ sin 2023 cos (^] F{<t>TA'i:0T:e'i)F{^i:.<t>'k;07:,0'f,), (43) 



B ^ sin2023sin (^^^ F(C, '^i^; 'A^; C ^D, (44) 

c ^ 4f(c, 01^; ^r, + S23^(c, 0?.; ^r, (45) 

Thus, we can calculate the coefficients A and B, which are related to the magnitude of the CP effect, by 
using our approximate formula. We can see the following from the expressions of A, B and C . The expres- 
sion of C is given as the sum of Fh and Fi, where we use the abbreviation Fh and Fi as i^(0™, 0^; 0™, 0)j) 
and F{(j)'^ 9^,9f). On the other hand, the expressions A and B are both given as the product of Fh 
and Fg and furthermore multiplied by the oscillating factor related to A32. This is the main difference 
between A, B and C. However, all the coefficients depend on the function F. In the following, we study 
the behavior of this function F. 
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At first, we divide F given in H88|l into three parts as 



i^(</)",0^;0'",0^) = F1+F2 + F3, (46) 

Fi = sin 20™ cos sin 20™, (47) 

F2 = cos^0"sin0'=sin2r, (48) 

F3 = sin2(/)'"sin(/)'=sin(26''=-46i"). (49) 

This separation of the function F is useful to understand, which contribution becomes large in the ampli- 
tude because Fi, F2 and F3 correspond to the MSW effect in the mantle, and in the core, and the mantle- 
core effect, respectively. By using the above expressions, the following interpretation in refs. jl2l can 
be understood more clearly. 

1. cos 20™ = and sin^'^ = 

Only Fi remains and the function takes the form F = ± sin 29™ because F2 — F^ = due to 
sin^'^ = 0. In the case that the above conditions are approximately satisfied around the MSW 
resonance region of the mantle, namely around the energy determined by sin 20™ = ±1, the function 
F is enhanced. 

2. sin0™ = and cos cj)'^ = 

Only F2 remains and the function takes the form F = ± sin 20^ because Fi = F3 = due to 
sin0™ = and cos(/)^ = 0. In the case that the above conditions are approximately satisfied around 
the MSW resonance region of the core, namely around the energy determined by sin 20^ = ±1, the 
function F is enhanced. 

3. COS0™ = and cos0^ — 

Only F3 remains and the function takes the form F = ± sin(20° — 40™) because Fi = F2 = due to 
cos 0™ = and cos (p'^ = 0. Around the energy determined by sin(20'^ — 40™) = ±1, the function F is 
enhanced. This can be large, even if both effective mixing angles in the mantle and in the core, 0™ 
and 0^, are small. It is considered as the mantle-core effect. It is realized in the case that A31 takes 
the intermediate value of the matter potentials a™ and a"^, respectively, for the mantle and the core. 

4.3 Interpretation of Numerical Results 

In this subsection, the numerical result for L = 12000km can be understood, by using the analytical 
expression derived in the previous subsection. All the coefficients A, B and C are determined by the 
functions Fg and Fh- Here, we study the behavior of Fi and Fh in the energy region larger than E — 1 
GeV. We can approximate Fi by using the fact A21 ^ a at F > 1 GeV. That is, the oscillation part and 
the mixing angle are approximated by 

IS.1L aL 

91 = - ~^ ^ const. (50) 

. A21 sin 2012 A21 sin 2012 1 
sm20^ = — ~ cx — (51) 

V21 cos 2012 - a)2 + All sin^ 20i2 ° ^ 

from (|15() . As a result, we can also approximate Fg from (|46|I - H49|I as 

Fi (X ^. (52) 
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Thus, the value of decreases proportional to the inverse of the neutrino energy. 

On the other hand, some of the peaks appear in Fh corresponding to Fhi, Fh2 and Fhz, since Fh 
includes the 1-3 MSW effect in the considered energy range. Fig. 5 shows the component of Fh by 
using our analytical expression (|46|l . where we use the matter densities in the mantle and the core as 
p-m _ 4,7g/cin'^ and p'^ = ll.Og/cm'^, the electron fraction as yj" = 0.494 and = 0.466, calculated by 
the PREM in the case of the baseline L = 12000 km. 



Component of F^ (sin^ 26*13 = 0.10) 



Component of F^ (sin^ 26*13 = 0.04) 
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Fig. 5 Components of F^ calculated from our analytical formula with sin^ 29i3 = 0.10 and 0.04 from left 
to right. Solid, dashed, dash-dotted and dotted lines correspond to F^, F^^, F^^ and F^^, respectively. 

In these figures, the solid line shows the magnitude of F^, and the dashed, dash-dotted and dotted 
lines show the magnitude of i^^i, ^^'^ ^h3^ respectively. Furthermore, we represent the values of each 
component F^i, Fh2 and Fh^ at three peaks in Table 2. 



sin^ 26*13 


E (GeV) 


Fh 


Fhi 


Fh2 


Fhs 


0.10 


2.1 


0.821 


0.278 


0.509 


0.034 


0.10 


3.0 


0.726 


-0.029 


0.011 


0.744 


0.10 


5.4 


-0.810 


-0.821 


0.015 


-0.004 


0.04 


2.0 


0.648 


0.093 


0.532 


0.023 


0.04 


3.1 


0.999 


0.051 


0.055 


0.893 


0.04 


5.7 


-0.535 


-0.545 


0.013 


-0.003 



Table 2. Components of Fh calculated from our analytical formula with sin^ 2(?i3 = 0.10, and 0.04. 

Fig. 5 and Table 2 show that the peak in the right hand side is dominated by Fhi and mainly depends 
on the MSW effect in the mantle. The MSW resonance in the mantle is realized at the condition a™ = 
A31 cos 26*13. The energy determined by this condition \s E ^ '^2'./lGjV"'' ^ ^''^ P^^^ 
hand side is dominated by Fh2 and mainly depends on the MSW effect in the core. The MSW resonance 
in the core is realized at the condition a"^ = A31 cos 26*13. Noticing the relation ~ 2.5 x a™, the peak 
energy is given by around E ^ 5.7/2.5 ~ 2.3 GeV. The energy of these peaks do not largely depend on 
the value ^13 in the case of sin^ 20i3 <C 1. Furthermore, it is shown that the mantle-core effect mainly 
contributes to the peak at the center, when Fh^ becomes large. The energy determined by the condition 
sin(26*'= - 46*") ~ 1 is about E = 3-4 GeV for sin^ 26*13 = 0.04. In the case of sin^ 26*13 = 0.10, this 
condition cannot be satisfied in any energy region and as a result the enhancement is weakened. This 
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phenomena is interesting because the value of for sin^ 2013 = 0.04 (small mixing) is larger than that for 
sin^ 2^13 = 0.10 (large mixing). It is interpreted as the total neutrino conversion pointed out by Petcov et 
al. [12 . 

Next, let us study how we can understand the behavior of A, B and C. ^From 10^, C is approximated 

by 

C=\{Ff+Fl)^\Fl (53) 

where we neglect , because of its smallness compared to F^ as shown in fig. 5. Actually, the C- function 
has almost half of the size of the F^-function. Therefore, C has three peaks as F^. P{ve v^) in ref. 
[T^ corresponds to C in this paper. It means that the terms A and B, related to the CP phase, were not 
considered in previous papers. 

Next, we obtain the expressions for A and B from (|47|l and 148|l as 

A ~ cos [ F^F,, cx -cos f — Fft, (54) 

B ~ sm ( — ^ j FeFh oc — sm ( — ^ 1 F^ . (55) 

^From these expressions, we can see the following. First, the mantle-core effect, which is different from the 
usual MSW effect, appears not only in C but also in A and B because of the multiplication of F/j. Second, 
A and B are suppressed compared with C in the energy range F > 1 GeV because of F^ ex 1/F. Third, 
A and B depend on the oscillation part -^^^ additionally to LjE dependence included in F. Because of 
this factor, the oscillation phases of A and B have a difference of about a quarter of the wavelength. 

5 Summary 

In this paper, we investigate the matter effect included in the terms related to the CP phase, particularly 
in the case that neutrinos pass through the earth core. The results are summarized as follows. 

1. Our approximate formulas ©-Q include non-perturbative effect of the small parameters a = 
A21/A31 and S13. As a result, the precision of the formula is rather improved compared to the 
formula which includes up to second order of a and S13 around the MSW resonance regions. 

2. We numerically calculate the coefficients A, B and C for the baseline length L = 6000 km and 12000 
km by using the PREM as the earth matter density. As a result, the magnitude of A and B can reach 
a few ten % of C around the three main peaks of C even for the case of including non-perturbative 
effect. 

3. We give the qualitative understanding of the behavior for A, B and C by using our approximate 
formula. The mantle-core effect, which is different from the usual MSW effect, appears not only in 
C but also in A and B, although the effect is weakened. 

^From the results of this paper, it has been found that the effects of the leptonic CP phase can be 
comparatively large in the oscillation probability, when neutrinos pass through the earth. We should 
consider the CP phase effects more seriously in order to extract the information on ^13 and the sign of 
Amgj^ in future experiments. 
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